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Abstract 

We calculate the exact analytical coefficients of the /3-expansion of the grand canon- 
ical partition function of the unidimensional Hubbard model up to order /3 5 , using an al- 
ternative method, based on properties of the Grassmann algebra. The results derived are 
non-perturbative and no restrictions on the set of parameters that characterize the model are 
^ required. By applying this method we obtain analytical results for the themodynamical quan- 
tities, in the high-temperature limits, for arbitrary density of electrons in the unidimensional 
chain. 
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1. Introduction 

Working with fermionic variables seems discouraging, for their non-commuting nature, 
in a certain way, compells one to a higher degree of care than that required by commutative 
variables. Nevertheless, Grassmann algebra properties justify their use in many circumstances 
[1]. With that in mind, we have recently described an alternative method to calculate the 
terms of the j3- expansion of the grand canonical partition function of periodic unidimensional 
self-interacting fermionic models, in which Grassmann algebra properties play a central role. 
No auxiliary fields are needed, and we express our results in terms of matrices with com- 
muting elements. The general approach for a periodic unidimensional fermionic model has 
been applied to the unidimensional Hubbard model up to order (3 s [2]. An important point 
about the method developed in reference [2] is that even though the unidimensional Hub- 
bard model has exact solutions [3], the analytical expressions are only known in the half- 
filled case. This drawback hinders the analytical evaluation of the partition function for the 
model from the knowledge of its energy spectrum when we are not considering the half-filling 
case. Takahashi[4] derived a closed expression for the grand canonical partition function of 
the unidimensional Hubbard model, but besides requiring the formulation of some additional 
hypotesis, a simple closed expression for the grand canonical partition function is only ob- 
tained in the strong coupling limit. The literature offers many examples of high temperature 
expansions of the grand canonical partition function for the Hubbard model in different space 
dimensions, some of them up to order (3 9 [5]. However, all these works referer to either some 
approximation in which one of the characteristic constants of the model has to be much bigger 
than the other, or some consideration based on important numerical analysis. 

Our results are analytical and do not rest upon any additional hypothesis on the constants 
that characterize the model. We should point out that we do not perform a perturbative 
expansion valid in the high temperature limit [5], but a /3-expansion of the grand canonical 
partition function [6] where we calculate the exact analytical coefficients of the terms up to 
order (3 5 for any density of electrons in the unidimensional chain. 

In the present paper, we will apply the approach developed in reference [2] to get the exact 
terms at orders (3 A and (3 5 of the grand canonical partition function of the unidimensional 
Hubbard model. In section 2 we present a review of the results of reference [2]. In section 3, 
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we write down the unidimensional Hubbard model and the Grassmann functions necessary in 
the calculations that follow. In section 4 we obtain the coefficients at orders (3 4 and (3 5 of the 
/3-expansion of the grand canonical partition function of the unidimensional Hubbard model. 
We use the grand potencial derived up to order /3 4 to calculate some physical quantities. A 
certain property of the multivariable Grassmann integrals — namely, its factorization — opens 
the way to extending the method of reference [2] to orders higher than (3 3 . This factorization 
property is presented in Appendix A through one example. In Appendix B we introduce a 
graphical notation that greatly simplifies the calculations. In Appendix C, we present a table 
of the necessary multivariable Grassmann integrals. Section 5 is dedicated to our conclusions. 



2. Review of Previous Results [2] 

The grand canonical partition function of any system is given by: 

Z(/?; M ) = Tr(e^ K ), (2.1) 
where (3 = k is the Boltzmann constant, T is the absolute temperature, and 

K = H-^/N, (2.1a) 

H is the hamiltonian of the system, [i is the chemical potential and N is the total number of 
particles operator. 

In the high temperature limit, (3 <C 1, Z((3; fx) has the expansion 

Z(P^) = Tt[1-PK}+ jrtPL Tr[K"], (2.2) 

n=2 

that we call the /3-expansion of the grand canonical partition function. 

For any self-interacting fermionic quantum system, in reference [2] we used the Grass- 
mann algebra to show that 

. 2nN 77j Au r)j 

Tr[K n ] = / J! dr]ldf]I eI ' J=1 X 
i=i 

x JC®(f}, rj;u = 0) K,®(fj, rr,u=l) ■■■ JC®(fj, rj;u = n-l), (2.3) 



where the matrix A is given by, 



A TT © 



(D A u 



(2.3a) 



and 



A.TT = A^ = 



(2.36) 



VljVxiV ©VxAf ©VxJV ••• liVxiV/ 

Each matrix A* 717 has dimension niV x nN, 1at x at been the identity matrix in dimension 
N x N and QvxTV the null matrix in this dimension. N is the number of space sites and n is 
the power of the (3 term. The non-null elements of A acr , a =f and a =|, are 



a n = 1, 



7 = l,2,---,niV 



yu , a 7 , 7+iV = -l, 7=1,2,- ..,(n-2)JV 



4 OCT _ 



(2.3c) 



I fl(n-l)iV+/,7 



1, 7=1,2, 



AT. 



The Grassmann function )C®(fj, rj) is the kernel of the fermionic operator K in the normal 
order [7]. In writing down eq.(2.3), we have used a particular mapping for the Grassmann 
generators [2], that greatly simplify our calculations: 



and 



7^{xi,T v ) = IJvN+l 



(2.4a) 



(2.46) 



where I = 1, 2, • • • , N, and v = 0, 1, • • • , n — 1. The mappings (2.4a-b) can be summarized as 
[8]: 



Va{xi,T u ) = V[il - zln+u]N+l 

The generators f) a {xi,T v ) have an equivalent mapping. 
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(2.4c) 



In reference [9] we showed that the Grassmann integrals (2.3) can be written as co-factors 
of the matrix A' 70 ". As we have stated before, the calculation of the co-factors of matrix A " " 
gets simpler when we diagonalize it through a similarity transformation 



where the matrix D is, 



D = 



p-l Af rap = D 



I AiIatxAT Qvxtv 



QvxTV 



(2.5) 



(2.5a) 



Xi, i = 1,2, •••,n, are the eigenvalues of matrices A aa , a =T>i- At the same time, we 
transform of the anti-commuting variables, 



7/ = P 1 r) and r}' = r)~P. (2.56) 

The matrices P and P _1 have a block-structure, and in reference [2] we got the elements 
of matrix P and its inverse for the particular cases n = 2 and n = 3, besides the eigenvalues 
of matrices A " ", a =T, !• However, for any value of n, we have that 



and 



with is, v' = 0, 1, 



and 



P 



(n) 



(») = J =e ^(2^' + 1)(^+1)_ 



n 



f , £(2i/' + 1)(i/+1) 



, n — 1, and 



\P n ^l j0 liVxJV ••• Pn-l^-l^NxN / 

/ ( n )n ( n ) n \ 

' Qqo %-NxN ■■■ ^O.n-l'WxfV \ 

L (n) -|, (n) -|, , 

\g n _ 1)0 ljVxJV ••• ? n _i in _ilArxAf / 



(2.6a) 



(2.66) 



(2.6c) 



(2.6a 1 ) 
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The diagonal elements of matrix D are: 

A W = i_ e ^+i) ) i/ = 0,l,---,ra-l. (2.6e) 

where the eigenvalues are iV-fold degenerated, N being the number of space sites. Due to 

(n) (n) 

space translation symmetry, we should note that the elements p u J, and q u J, do not carry any 
space site index. 

Once the matrix A has a block-structure as depicted in eq.(2.3), the integrals (2.3) are 
equal to the product of the integral in the sector a a =|T times the integral in the sector 
oo =H . The Grassmann integrals in the sector || have the form: 

nN 

/nN ^2 fjiAlJrij 

Y[dr H dfj i f}i 1 r ]kl ---fii m i lkm , (2.7) 

i=l 

with L = {Zi, • • • , l m } and K = {k\, • • • , k m }. The products 7777 are ordered in such a way 
that h < h < • • • < l m and fci < k 2 < ■ ■ ■ < k m . From reference [9], the result of this type of 
integrals is equal to: 

M(L, K) = (_i)('i+i2+-+'™)+(fci+fc 2 +-+fcm)^( ij K ), (2.7a) 

where A(L, K) is the determinant of the matrix obtained from matrix A by deleting the lines 
{hi ■ • • 5 ^m} and the columns: {ki, ■ ■ ■ , fc m }. The Grassmann integrals to be calculated in 
sector || are of the same type as that of eq.(2.7). 

After the similarity transformation (2.5) and the change of variables (2.5b), in a schematic 
way, the integral (2.7) becomes 



nN 

/nN ViDijVj 
i=l 

(2.8) 

where D is a diagonal matrix whose entries are given by eq.(2.6e). 

We should point out that the relations (2.6a)-(2.6e) are valid for any unidimensional 
self-interacting fermionic model with space translation symmetry. 
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We have a large number of integrals that contribute to eq.(2.3). For a discussion on some 
useful symmetries and their use in the reduction of the number of contributing integrals, the 
reader is referred to [2]. 

3. Unidimensional Hubbard Model 

The hamiltonian that describes the Hubbard model in one space dimension is [10]: 

N N N 

» = 1 <T = -1,1 1=1 i=l <T = -1,1 

where aj^ is the creation operator of an electron in site % with spin a and a^ is the destruction 
operator of an electron in site i with spin a. The first term on the r.h.s. of eq.(3.1) is the 
kinetic energy operator. All diagonal elements of Uj are equal, tu = Eo, the only non-null 
off-diagonal terms are U^-i = U^+i = t, where i = 1, 2, . . . , N, and they contribute to the 
hopping term. U is the strength of the interaction between the electrons in the same site but 
with different spins. We have defined A^ = — ^gusB, where g is the Lande's factor, \xb is 
the Bohr's magneton and B is the constant external magnetic field in the z direction. 

The periodic boundary condition in space is implemented by imposing that ao<7 = &-n<t 
and a7v + i jCr = a.\ a . Therefore, the hopping terms tioaj^aoo- and tA^jv+iaj^ajv+i,^ become 
tiN&itj&Na and tjv,i a jvcr a icr respectively. We point out that the hamiltonian (3.1) is already 
in normal order. 

The kernel of the operator K (eq.(2.1a)) for the unidimensional Hubbard model on a 
lattice with N space sites, written in terms of the generators fji and rjj, is equal to 

N 

3C®(f], 77; U)) = ( E ° + aXB - V { (l^l n+l/]N+l \^l n+v]N+l + 

1 = 1 cr=±l 

N 

1 = 1 cr=±l 
N 

+ ^ UV(n+u)N+l V(n,+v)N+l VvN+l TfuN+l, (3.2a) 
1=1 

with the periodic spatial boundary condition: 
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t V[(l^l n+U ] N+N V[(lz-zl n+J/ ] N+N+1 - t Vi±_zl n+iy]N+N V[il_H n+iy]N+1 (3-26) 

and 

t V { ii^l n+v]N+1 V[(l^l n+iy]N = t V { (l_zl n+iy]N+1 V { (l^l n+u]N+N , (3.2c) 
and the anti-periodic boundary condition in v: 

V [ (i^i n+n]N+ i = -V { ii^l n]N+l , (3.2d) 

for I = 1,2, • • • , AT, and a =T>!- We have used the mapping (2.4c) to write the previous 
expressions. 

In order to write down the terms that contribute to Tr[K 4 ] and Tr[K 5 ] in a simplified 
way , we define: 



JV 

£(fj, 77; i^;a) = J2 V { (i^l n+I/]N+l V { (i^l n+I/]N+l 5 (3.3a) 
1=1 

N 

T- (77, rj; v;a) = J2 V[i^l n+v]N+l V [ ii^l n+v]N+l+1 5 (3-36) 



and 



2 

1=1 



N 



T+(fj, ry; u; a) = ^ ^a^)^^ ^(w)^^^ . (3.3c) 
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We also define 



£(ri,rr,u)= E (v)£(v,V; (3.4a) 

<7 = ±1 

T~(v,V,v) = tT~(fj,r];u;a), (3.46) 

<7=±1 

T + (f},r];v) = Y tr+ %m ^; (7 ) ) (3-4c) 



=±1 



and 



where E(a) = E — a\s — The term £(fj, 77; v) represents the diagonal part of the kinetic 
energy , T~ (77, 77; v) and T + (fj, r\\v) are the hopping terms and U(fj,r];i') is the fermionic 
interaction term. 

For the unidimensional Hubbard model, the Grassmannian function JC®(fj, 77; v) is written 

as 

K®%7]\v) = £(77,77;*/) +T-(fj,ri;v) + T + (f},ri;is) + U(fj,ri;u). (3.5) 

4. The Exact Coefficients of the (3— Expansion of the Grand Canonical Partition 
Function for the Unidimensional Hubbard Model 

In eq.(2.2) we have the /3-expansion of the grand canonical partition function for any 
quantum system. For the particular case of self-interacting unidimensional fermionic models, 
the expression of Tr[K n ] is given by eq.(2.3). 

In reference [2] , we calculated the exact coefficients of the terms (3 2 and /3 s of the expres- 
sion (2.2) for the unidimensional Hubbard model for arbitrary values of the constants Eq, t, 
U and (j,, that characterize the model, and for any value of the constant external magnetic 
field B. 

The evaluation of integrals has been performed by a number of procedures (computer 
programs) developed by the authors in the symbolic language Maple V.3, that consist in the 
computational implementation of the method described in reference [2]. We have called this 
package of procedures GINT. 

The method applied in [2] greatly simplifies the calculations made in reference [11], 
but memory utilization problems have appeared as we tried to go beyond n > 3. Luckily, 
the factorization property of multivariable integrals of type (2.8) allowed us to optimize the 
performance of the package. In Appendix A we consider one typical Grassmann integral to 
exemplify the factorization of the sub-graphs. 

The procedure perm is one of the procedures contained in the package, being a useful 
tool to calculate the independent non-null terms [12] that contribute to Tr[K 4 ] and Tr[K 5 ]. 
In this procedure are implemented the symmetries discussed in reference [2]. The procedure 
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gint has been used to calculate the multivariable Grassmann integrals, taking into account 
the factorization into sub-graphs. The package can be downloaded through ftp from the site 
http :/ www. if. uff. br. 



4.1. Calculation of Tr[K 4 ] 

For the case n = 4, we get from eq.(2.3) that 

Tr[K 4 ] = / J"[ e 7 - J=1 x 

1=1 

x K®% 77; v = 0) K®(fj, 77; v = 1) /C®(77, t/; 1/ = 2) K®% 77; 1/ = 3). (4.1.1) 

The expressions of A' 70 ", p^l, and are obtained from eqs.(2.3b), (2.6a), (2.6b) and 
(2.6e). 

From eq.(3.5), for the unidimensional Hubbard model, we have that the Grassmann 
function K,®(f), 77; v) is equal to 

K®% 77; v) = S(fj,r);v) + T~(fj,r];u) + T + (fj,rj;u) +U(fj,r};u). (4.1.1a) 
We defined a simplified notation, 



/2nN ^ 77/ Au rjj 

1=1 

x 0i fo77;z/i)---<9 m (77,77; z/ m ) (4.1.2a) 



and 



/(3-o-)niV ^ »7(i-<7)nJV+i Arj '7 ( i_ ct) „ jv+ j 

JJ e J - 7=1 x 

I=(l-a)nN+l 

x 01(77,77; 1/1) •••<9 m (77,77;^), (4.1.26) 

where Oj(fj,r];vj) are Grassmann functions. 

The independent terms that contribute to Tr[K 4 ] are: 
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Tr[K 4 ] =< £ ,£o,£o,£o > +4 < U, £ , £q, £q > +2 < U,£ M,£o > +8 < W,T~,T + ,£ > + 
+ 4 < U,U,£ ,£o > +A<UMM,£o > + <UMMM > +4 < T-,£ ,T+,£ > + 
+ 8 < T-,W,T+,£ > +4 < T~,U,T + ,U > +8 < T-,T+,£ ,£o > + 
+ 8 < T-,T+,W,£ > +8 < T~,T+,U,U > +2 < T~ ,T + ,T~ ,T + > + 
+ 4 < T-,T-,T + ,T+ > . (4.1.3) 

In order to calculate the terms on the r.h.s. of (4.1.3), we need the result of a set of 
Grassmann multivariable integrals that are presented on Appendix C. We used the procedure 
gint to calculate them. Before using the procedure gint to calculate the terms in eq.(4.1.3) 
that include £q, T~ and T + , we need to explicit the contributions coming from the sectors 
a =| and a =|. For example, in the term < £ (0), £ (1), £(2), £ (3) > we have 16 terms when 
we explicitly write down the spin-sectors. However, the number of terms is diminished when 
we use the symmetries discussed in reference [2] and the fact that = The application 
of these symmetries is simplified by using the graphic notation explained in Appendix B. 

By taking into account the results of integrals in Appendix C and their contributions to 
the sum over the space indices in each term of Tr[K 4 ] (eq.(4.1.3)), we finally obtain 

Tr[K 4 ] = N 4 (UAE 3 + — + -U 2 AE 2 + AE A + —U 3 AE) + N 3 [3AE 2 \ 2 B + 3AE 4 + 
\ 256 8 16 /V 

+ 3t 2 UAE + -UAE 3 + -U 2 X 2 B + 6t 2 AE 2 + — U 4 + ^U 2 AE 2 + 
2 16 128 16 

+ ^t 2 U 2 + -U 3 AE + -UAEX 2 B ) +N 2 ( -AE 4 + —U 2 X 2 B + 3UAE 3 + —U 2 AE 2 + 
8 4 2 B ) \A 16 B 16 

+ 3t 2 UAE + \aE 2 X 2 b + 3t 2 X 2 B + h 4 B + lt 2 U 2 + M 4 + M 2 AE 2 + ^U 3 AE + ^-U 4 ] + 
2 4 8 16 256 

-n( 3t 2 X 2 B + -t 4 + 3t 2 AE 2 + -UAE 3 + -t 2 U 2 + -U 2 X 2 B + - A AE 4 + 
\ 2 2 4 8 4 

+ 3t 2 UAE + \aE 2 X 2 b + \u 2 AE 2 + -X 4 B + \u 3 AE + \uAEX 2 B + ^-U 4 
2 8 4 8 2 128 

(4.1.4) 

We use the short notation: AE = Eq — /x. 
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4.2. Calculation of Tr[K 5 ] 

For n = 5, the expression of Tr[K 5 ] obtained from the procedure perm, where each term 
is multiplied by the respective constant, is: 

Tr[K 5 ] = 5 <U, S , £ , So, So > +5 < U, S , U, S , S > +5 < U, U, S , S , S > + 
+ 5 <UM,SoM,So > +10 < U,U,T-,T + ,So > +5 <UMM,So,S > + 
+ 5 <UMMM,So > + <UMMMM > + < S ,So,So,S ,So > + 
+ 10 < T-,T+,T-,T+,£ > +10 < T-,T + ,T-,T + ,l( > +10 < T~MM,T + ,So > + 
+ 10 < T-,T + ,So,U,So > +10 < T-,T+ So,S ,S > +10 < T", T+, W, £ , £ > + 
+ 10 < T-,T + ,U,U,S > +10 < T~,T + ,U,U,U > +10 < T-,S ,T + ,U,So > + 
+ 10 < T-,So,T+,S ,So > +10 < T-,U,T + ,So,S > +10 < T~ ,U,T+ ,U,£ > + 
+ 10 < T~M,T + MM > +10 < T + ,r-,T-,T+,^o > +10 < T + ,T-,T-,T+M > + 
+ 10 <U,T- M,T + ,So > +10 < W,T-,£ ,T+,£ > +10 < U,T~ ,T + ,S ,S > + 
+ 10 < U,T- ,T + M,So > +10 < T-,T-,T + ,T+,£ > +10 < T~ ,T~ ,T + ,T + ,U >, 

(4.2.1) 

where we are using the convention (4.1.2a) to write down each term. 

In n = 5 we have seven new type of integrals that do not have an equivalent one for 
n < 5. We table those integrals in section C.2 of Appendix C. 

After a long but convergent calculation, we get 

Tr[K 5 ] = N 5 {— A U b + — U 4 AE + -UAE A + — U 3 AE 2 + AE 5 + -^AE 3 ] + 
y 1024 256 4 32 8 J 

+ N 4 ^u 2 AE 3 + — U 2 AEX 2 B + — C/ 3 t 2 + 10t 2 AE 3 + 5A£ 3 A 2 + 5A£ 5 + — UAE A + 
V 16 16 B 32 s 4 

+ —UAE 2 X% + -t^ 2 A£ 2 + — U 5 + -f/ 3 A 2 + — U 3 AE 2 + -^-£/ 4 A£+ 
4 B 2 512 64 B 64 128 

+ -U 2 t 2 AE] + N 3 ( 15t 2 AE 3 + 15t 2 AEX% + — £/A£ 4 + -UAE 2 X 2 B + 
8 J \ 16 8 

+ l ^UX% + l ^AE 3 X 2 B + ^AE 5 + ^A£A 4 B + ±jW + ^[/ 2 A£ 3 + ^[/ 2 A£;A|+ 

+ ^Ut 2 X 2 B + 7 4ut 2 AE 2 + 15t 4 A£ + ^U 4 AE + JVa 2 ? + ^-U 3 AE 2 + 
4 ' 4 256 32 2 

+ -U 2 t 2 AE + — £/ 5 + -[/ 3 t 2 ^| + N 2 f-^U 2 t 2 AE + -£/ 2 A£ 3 - 
8 1024 32 / V 4 4 
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^-U 2 AEX 2 B + ^-U 5 - ^UAE* - 7 4uAE 2 X 2 B - ^UX% - ^t 4 AE- 



512 16 8 " a 16 a 2 

^rfAE 2 - -Ut 2 X 2 B + ^U*AE--AE 3 X 2 B - -AE 5 - -AEX% - \ ~>r.M:/ 
2 5 128 2 B 4 I 

- 15t 2 AEX 2 B - f A U*X 2 B + ^U'AE 2 - fut 4 - \uh 2 ] + 

+ n(--U 2 AE 3 - -Ut 2 AE 2 + -Ut 2 X 2 B - -U 2 t 2 AE - -U 4 AE - -U 3 AE 2 + 
\ 2 2 2 ^ 2 32 32 

+ lt/ 3 A| - — U 5 - -UH 2 - -UAE^ + -UX% 
32 B 128 8 4 4 B 



(4.2.2) 



We continue to use the notation: AE = E — /j,. 

Certainly, the most subtle part in calculating expression (4.2.2) comes from the product 
of Grassmann integrals for different a-sectors when we have Grassmann generators at the 
same space indice. In this case, we have to suit the conditions satisfied by both integrals and 
calculate the contribution of the product to the sum over space indices. 



4.3. The (3— Expansion of the Grand Potential Up to Order /5 4 and 
Physical Quantities 

The relation between the grand potential W(/3; //) and the grand canonical partition 
function Z(/3; fx) is 

W((3;tx) = ~lnZ(P;Li). (4.3.1) 
The /3-expansion of Z{J5\ /x) up to order f3 5 is (see eq.(2.2)): 

Z((3,li) « Tr[l-/3K] + ^Tr[K 2 ] - ^Tr[K 3 ] + ^Tr[K 4 ] - ^Tr[K 5 ]. (4.3.2) 

The first term on the r.h.s. of (4.3.2) was calculated in reference [11], the second and third 
terms were calculated in reference [2] and in its last two terms we substitute the results of 
eqs. (4.1.4) and (4.2.2). 

From eqs. (4. 3.1) and (4.3.2), we get the grand potential up to order /3 4 , that is, 
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768 



/*) = -JV{| In 2 + ("^^ 2 A| + ±-U 2 t 2 AE + ^U 5 + ^Ut 2 AE 

_ — U 3 X 2 - —UX% + — UAE 4 + — U 4 AE + —U 2 AE 3 + —U 3 t 2 )/3 4 - 
768 B 96 B 96 768 48 64 J f 

- ( -t 2 UAE + — UAEX 2 + —AE 4 + —t 4 + —X 4 + —U 4 + —UAE 3 + 
\8 16 B 96 16 96 B 1024 48 

+ -t 2 A£ 2 + — U 3 AE + —t 2 U 2 + — U 2 X 2 B + —U 2 AE 2 + —AE 2 X 2 R + 
8 192 96 64 B 64 16 5 

+ -t 2 X 2 B )p 3 + (-— + -UX% - —AE 2 U - -AEU 2 )/3 2 + 
8 B ) H V 64 16 16 16 y 

+ (I AE 2 + \^B + l^U+^ + ^U 2 y- (AE +^)+ O(^) } . 

(4.3.3) 

It is important to stress out that the coefficients of the /3-expansion of function W(/3; n) are 
exact for any set of constants: Eq, t, U, [i and B for the unidimensional Hubbard model. 
From expression (4.3.3) we can get the strong limit approximation by taking U 3> t, as well 
the atomic limit approximation when U < t. 

From expression (4.3.3), we can derive any physical quantity for the model at thermal 
equilibrium at high temperature. As examples, we consider the following quantities: 

i) specific heat at constant length and constant number of fermions: Cl(/3). 



d/3 

where k is the Boltzmann constant. From eq. (4.3.3), we get, 



(4.3.4) 



C L ([3) = Nkl ( —U 5 - —UX% + -UAE 4 + — U 3 AE 2 + -U 2 AE 3 - —U 3 X 2 B - 
V 1 \V256 24 B 24 192 12 192 B 

_ ^ut 2 X 2 + -Ut 2 AE 2 + -U 2 t 2 AE + — U 4 AE + -U 3 t 2 ) [3 5 + (-—U 4 - 
4 B 4 4 192 16 / V 256 

- \t 2 AE 2 - \X 4 B - \t 4 - \AE 4 - \t 2 U 2 - ^U 3 AE - \uAEX 2 B - \aE 2 X 2 b - 

~ h 2 X 2 B - lt 2 UAE-^-U 2 X 2 B - ^U 2 AE 2 - -UAE 3 ](3 4 + 
11 16 16 4 J 
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+ (-lu 2 AE - ^U 3 + S -UXl - \lJAE*y + 

+ (hjAE + t 2 + ^U 2 + 1a| + \AE 2 \ ii 2 + G(/3 6 ) [>. (4.3.4a) 



ii) average energy per site: < h > (/?). 

The simplest way to derive the average energy per site from the grand potential, is to 
scale the constants: (E , t, U,Xb) — > (ctE , at, all, olXb) and substitute in eq.(3.2a) to obtain 
W(/3; //; a). From the scaled grand potential, we have 



<h> ian|M 



From eq. (4.3.3), we get that 



< > 09) = E + j + (-t 2 - ^U 2 - l -E U + hj^ - \e 2 + ±E of i - \x 2 B y+ 
+ &oU 2 - \u 2 » - ^UXl + ^UE 2 - \UE , + ±U» 2 + A^, 2+ 

+ (- ^ + ^ + \* + + k** + y/ u2 + 

+ Y6 U2E2 + 3^ V + T2 UE ° -T8 U ^ + 4^° - + k E «~ 

- ^U 2 E 0f i + l -UE 0l i 2 + l -Ut 2 E - ^Ut 2 ^ - l -E^ + \E 2 ^ 2 - 

- ^E fi 3 + \t 2 E 2 + -t 2 ^ - ^t 2 E n + \UE X 2 B - h 2 B E n- 

+ (- ik u5 + -ws u * E ° ~ l uH2 + k u * + m*" ~ i u ^ + 
+ vk u3x * - k UE * - k u ^ - ws u ^ + k u ^ + \ UE ^~ 

+ ^U 3 E^ - ^UE 2 f + \u 2 E 2 ^ - ^U 2 E^ 2 + -^UEoV 3 - ^U 2 t 2 E + 
+ ^Ut 2 X 2 B - ^Ut 2 E 2 - ^Ut 2 ^ 2 + \t 2 U 2 ^ + \ut 2 E^-^U 3 A^+ 



+ 0(P 5 ). (4.3.5a) 



14 



Hi) difference between average numbers of spin up and spin down particles per site: 
< n-f > — < nj_ >. 

From the definition of the grand potential (eq.(4.3.1)), we have that 



< „, > («- < ni > (« = (4.3.6) 
Up to order /3 4 , we get from eq. (4.3.3) that, 

< rz T >(/?)-< n; >(/?) = { -(C//? + 4) + (A£ 2 + £/A£ + ±A| + ^ 2 + 2^/3 3 + 
+ f + Ut 2 + ^U 3 )p 4 + 0((3 5 )\. (4.3.6a) 



iv) average of the square of the magnetization per site: < m 2 > (/?). 



< w 2 > (/?) = ^§ < (n lT - n a ) 2 > 



,1 ,2 1 



d[i dU 

where B is the external magnetic field. From eq. (4.3.3), we obtain that 

< ml > (13) = jsVbU + \u/3 + (-\UE + - ±U 2 + \\% - ±E?,+ 



+ 



+ 



J_u± + ^u 2 e 2 + —u 2 ^ 2 + -ueI - -?V+ 

1536 384 384 P 24 24 P 



384 
1 



U 3 E Q 



384 ^ 32 



3 :C/ 2 A 2 _ 1 2 A 2 



(4.3.7) 



384 



1 



+ -UE^ 2 + ±Ut 2 E - \ut 2 ^ + ±E 4 + ^ - ±E$^ + ^E 2 » 2 - 

- y 2 e ^" + l f2E o + \*t ~ \ t2E ^ ~ W2 U2Eop )^ + W5) }' (4 - 3 - 7a) 



where g is the Lande's factor and \xb is the Bohr's magneton. 



v) magnetic susceptibility: %(/3). 
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xW = -(- 2 sm) jj—gg? 1 - (4-3-8) 



From eq. (4.3.3), we obtain that 



xW) = -{\gvB) 2 



ut 2 + 



-U\ 2 B + 



384 



U 3 /3 4 + 



+ h 2 B + 1 -uae+ 



+ i A £ 2 + -[/ 2 
8 32 



(4.3.8a) 



5. Conclusions 

With the implementation of the factorization into sub-graphs of the Grassmann multi- 
variable integrals, we can certainly go beyond the calculation of the term /3 4 of the /3-expansion 
of the grand potential of the unidimensional Hubbrad model. Even though the physics for 
U > and U < are different, the results of section 4.3 apply equally well for both cases. 
Recently, dos Santos and Thomaz have applied the results of reference [2] to calculate the 
/^-expansion of the grand canonical partition function of the extended unidimensional Hub- 
bard model up to orde /3 3 [16]. But the important point is that the present approach opens 
the possibility to calculate the first terms of the /3-expansion of the grand canonical partition 
function of the Hubbard model in two space dimensions, as well as of unidimensional models 
with impurities. We believe that improvements on the present approach will render a valuable 
tool for tackling with such problems. 
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Appendix A 
Factorization of Grassmannian Sub— Graphs 



For calculating the co-factors of matrices A " " , a =T, j, it helps to have the value of their 
determinant. For arbitrary n, the determinant of these matrices is equal to 



n-l 



detA c 



= [n*s 



n) 



N 



(A.1) 



u=0 



where N is the number of space sites and A^ are the TV-fold degenerated eigenvalue of A c 
From eq.(2.6e) we have that 



We should notice that A 
We define: 



(n) 

(n-l)-i/ 



A (n) =1 _ e -(2,+l)_ 

= A<< 



(A.2) 



n-l 



^ = II A - n) - 



i/=0 



For calculating we need to consider the cases n even and n odd separately. 
For n even, expression (A. 3) can be rewritten as, 



n-2 
2 



= n 



i/=0 



2 - 2cos(-(2i/ + 1)) 



n 



where the last equality is already known [13]. 

For n odd, expression (A. 3) can be rewritten as, 



n-l , x 

p(n) = 2 x 2 n-l J| sin j JL(2v+\)\. 



From reference [13], we have that [14]: 



n-l , v 

2n_i n sin u^ (2i/+i) =1 ' 

i/=0 ^ ' 



that substituted in eq. (A. 5) gives 



(A3) 



(A4) 



(A5) 



(A6) 
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p(n) = 2 

From the results (A. 4) and (A. 7), for any n, we get that 



(A.7) 



n-l 



detA c 



[ n*s 

v=Q 



n) 



N 



>iV 



(A.8) 



To present the factorization of the Grassmannian integrals, we consider an example and 
use the graphic notation explained in Appendix B. 

Let us consider the integral for fixed space indices l\ and Z3, that contributes to 

< 5 (T),£i(T),£2(T) >, 

47V 



, — T A Vj Vj 

l(h,h,h) = / II dVid-Vi e J ' J=1 Vh Vh VN+h VN+h V 2N +i 3 V2N+i 3 



i=i 



(A.9) 



where l\ 7^/3. Under the similarity transformation (2.5), eq.(A.9) becomes 



AN 



QviOPOti Iv^lPlr-^Vv^N +lx VriN+h Vv 2 N+h Vr 2 N+h 

«/ t 1 ,,~ — n 

4JV 

Vi D u V.i 



1=1 



T l .' r 2= 



X 



^3=0 



x ; > Qiy 3 2P2iy 3 vl 3 N+l 3 Vl 3 N+l 3 



2» £ 



\( 4 M 4 ) 

i^l,^2= ^1 ^2 

rj ,t 2 =0 



X 



6 



1/3=0 /v ^3 



A 



(4) 



(A10) 



where to write the second equality on the r.h.s. of eq.(A.lO), we used the result (A.8). 
By explicitly writting down the expressions, we see that, 



2" E 



v l >"2=0 

T\,T 2 =0 



x(4K(4) 



, 47V £ fjj AH V.I 

I Yl drudfj! e'- J=1 77^ 77^ 77^+^ VN+h- 

J T 1 



7=1 



(AlOa) 



and 
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4N 



E 



3 




[ [ drjidr)_ 



47V 



E a Vj vj 



U 3 =0 




7=1 



j=i 



(A.lOb) 
(All) 



Using the graphic representation of Appendix B, we write result (A. 10) as, 

The factorization (A. 11) comes directly from the fact that the matrix D is diagonal (see 
eq.(2.5a)) and the result (2.7a). Once the presence of Grassmann generators in the integrand 
of integrals (2.8) correspond to cutting lines and columns of matrix D,then only for cutts at 
the same space index and any ^-indices we get co-factors of matrix D that are non-zero. In 
summary, the factorization of the type (A. 11) always happens when two or more space indices 
are different. 

In a similar way and by the reasons discussed before, it is simple to show that X(Zi,^,-^), 
where all the space indices k, i = 1,2,3, are distinct, is easily written as: 
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Appendix B 

Graphic Notation of the Multivariable Grassmann Integrals 



To exemplify our graphic notation for the graphs that contribute to eq. (2.3), for fixed 
value of n, we consider some terms of Tr[K 4 ]. This graphic notation is very helpful when we 
apply the symmetries discussed in reference [2] to identify equivalent terms in Tr[K n ]. 

We present the graphic notation through examples and its application to the identification 
of equivalent integrals. 

1) 

<£(T,0),£(|,1),£(T,2),£(T,3) >=£(T) 3 £(|)x 

2) 

The constants i?(t) and E(l) are defined just below eq.(3.4d). In order to show that terms 
(B.l) and<(2($) 0)i&($qt^,(fee2^4Tt^ iftvafiH$^e©C|t)he integrals under a cyclic translation in 
the temperature parameter v in each sector a =f and a =|, separately. Therefore, 

<5(T,0)5a,l)5(T,2)^(T,3)>=<5a,0),5(T J l),5(T,2),f(T,3)> (5.3) 
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3) 

Due to the presence of the term U(0), the integrals in the two cr-sectors have one space index 

. . < W(0)T-(T,l)r+(T,2)5(T,3) >= Ut 2 E{])x 

I in common. 



4) 

The terms (B.4) and (B.5) are equal, up to a multiplicative factor, due to the fact that 

A n = A 6W(o)r-a,i)r^a,2)f a,3) >= u^e^x 

The graphic notation was used along all the calculations and permitted us to considerably 
reduce the number of terms that contribute to the expressions of Tr[K 4 ] and Tr[K 5 ]. 
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Appendix C 

Useful Multivariable Grassmann Integrals at n = 4 and n = 5 



C.l. Useful integrals for n=4 

We need the result of twelve integrals only, to calculate the terms that contribute to 
(4.1.3). In this Appendix, we present the value of these integrals according to the conditions 
satisfied by the space indices [15]. 

1) 

4JV 

, 47V £ fjr AH VJ 

j{ 4) (0 = J Y[d Vl dfji e'^ 1 f)im =2 N ~\ (C.l.l) 



1=1 



for 1 = 1,2, ---,N. 
2) 



. 4iV £ ^ AtJ ^ 

(Zi, Z 2 ) = / [[ drudru e 1 ^- 1 rj h rj h rj N+h rj N+h 

J 7 1 



1=1 



2 N ~\ h = h, Zi = l,2,---,iV 
2 N ~ 2 , h^l u h,h = l,2,---,N. 



3) 



, 47V £ ^ A}T VJ 

H (luh)= / 1J_ dr}idr}i e 7 - 7 - 1 Vh Vh+i Vn+i 2 Vn+i 2 -i 



1=1 



= 2 N - 2 ,l 2 = h + l, h = l,2,---,N. (C.1.3) 



4) 



(C.1.2) 



, 47V £ ^ A}T VJ 

H (h,h,h) = I \\dr}idqi e 1 -'- 1 Vh Vh Vn+i 2 Vn+i 2 +i r] 2N +i 3 V2N+i 3 -i 

J T — 1 



1=1 

)7V-3 

7V-2 i , (C.1.4) 



2"" d , Z 3 = Z 2 + 1 
2 Ar_2 , h = h = h — 1- 
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5) 



. 4iV £ ^ A]\ VJ 

H (h,h,h)= / [[ driidqi e 7 ^- 1 r/ Zl r/ Zl r] N+h r] N+h r] 2N +i 3 V2N+i 3 

J T 1 



6) 



7=1 



= {^ N n ~\ h = l2,OT,h = l 3 ,OT,l 2 = l 3 (C.1.5) 

2 , h = h = h- 



4JV 



(h>fe,fe) = I \\driidqi e 1 -'- 1 Vh Vh+i Vn+i 2 Vn+i 2 r] 2N +i 3 V2N+i 3 -i 



i=i 

)N-3 

N—2 i 7 „„j ; 7 i i 

(C.1.6) 



2"~\ h = h + l 
2 N ~ 2 , l 2 = h and l 3 = h + 1 



7) 



A (h,hih) = / II dr li d Vi e 7 " 7 - 1 r} h r} h r] N+h Vn+i 2 +i V2N+i 3 r]2N+i 3 -i 

J i=i 

2"" 2 , Z 1 = Z 3 = i 2 + l. ( } 



„ 47V ^ji 

l£\liM,hM) = / J] dr li df li e 7 '^ 1 

^ 7=1 



% % 777V+Z 2 VN+l 2 V2N+l 3 V2N+l 3 V3N+l 4 





r 2 iV ~ 4 


h 






h 


= < 


2 N ~ 2 , 


h 






h 




,2 N ~\ 


h 



(C.1.8) 



h = or ? au permutations with 3 equal space indices 
h = h = U- 
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9) 



1iV 

1 9 (h,h,h,k)= [I drjidrji e 1 -'- 1 rj h rji 1+1 r] N+ i 2 r} N +h-i V2N+i 3 V2N+i 3 +iX 

J 1=1 

X V3N+l 4 V3N+l 4 -l 



2 7V-4 
2 N-2 



, h = h + 1 and Z 4 = h + 1, or, Z 2 = Z 3 + 1 and £4 = h + 1 1 9) 
, h = h + 1 and Zi = Z3 and Z4 = Z 3 + 1. 



10) 



4t j\ 

^io (I1J2J3J4) = / || dry/dr// e 7 - 7 - 1 Vh Vh-i Vn+i 2 Vn+i 2 -i V2N+i 3 V2N+i 3 +i x 

^ 7=1 

X V3N+l 4 V3N+U + 1 



2 N ~ 4 , Z 2 = Z 3 + 1 and Zi = Z 4 + 1, or, Z 4 = Z 2 - 1 and Zi = Z 3 + 1 
= <( 2^, / 2 = Z 3 + 1 and h = Z 3 + 2 and Z 4 = k + 1 (C.1.10) 
Z2 = Z3 + 1 and /1 = Z3 and Z 4 = Z3 — 1. 



11) 



^11 (Zi,Z 2 ,Z 3 ,Z 4 ) = / || c%c% e 7 " 7 - 1 ^1^1 VN+hVN+h V2N+i 3 rj2N+i 3 +i x 

^ 7=1 

X fj 3N+ i 4 r/3AT + i 4 _i 



2^" 4 , Z 4 = Z 3 + 1 

2 Ar ~ 3 , / 2 = Z 3 + 1 and Z 4 = Z 3 + 1, or, Zi = Z 2 and Z 4 = Z 3 + 1, ,^ ^ ^ 

or, li = l 3 + 1 and Z 4 = / 3 + 1 
2 JNr-2 J Zi = Z 3 + 1 and Z 2 = Z 3 + 1 and Z 4 = Z 3 + 1. 



12) 



,47V £ ^AtT^ 

^12 (^i,^2,Z 3 ,Z4) = / II rfy 7/d?/7 e 7 ^- 1 77^ 77jv+/ 2 ??tv+z 2 +i ??27V+z 3 V2N+i 3 x 

7=1 

X ^37V+Z 4 ^iV+Lt-l 



2^- 4 , Z 4 = Z 2 + 1 

2 7V_3 , Zi = Z 3 and Z 4 = Z 2 + 1, or, Zi = Z 2 + 1 and Z 4 = Z 2 + 1, 

or, Z 2 = Z3 and Z 4 = Z 2 + 1 
2 N ~ 2 , Zi = Z 2 + 1 and Z 3 = Z 2 and l 4 



(C.1.12) 



Z2 + I. 
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C.2. Useful integrals for n=5 

We present here the seven integrals that have no equivalent ones for n < 5; i.e., integrals 
that cannot be factorized into any of the integrals for n < 5 for all the conditions satisfied by 
the apace indices. In some graphs, we do not have generators in the integrand of integrals of 
type (2.7) at a given value of u, as we can see for example in the graphs presented in Appendix 
B. Those rings in the integrals of type (2.7) that have no associated Grassmann generators 
in the integrand, we call empty rings. For example, in (B.l) we have one empty ring at a =| 
{y = 1), and three empty rings at a =j {y = 0, 1, and 3). The integrals for n = 5 with empty 
rings give the same results to the equivalent integrals for n = 4. We have not demonstrated 
this property in general form for any n, but we have detected it by evaluating these integrals 
through the procedure gint. 

The seven integrals for n = 5 and the conditions satisfied by the space indices [15] are: 

1) 

5JV 

Q[ (h,h,h,l4,h) = jWdrndrn e 1 ^- 1 Vh Vh Vn+i 2 Vn+i 2 V2N+i 3 V2N+i 3 x 

x V3N+U V3N+U V4N+l 5 V4:N+l 5 

(2 N ~ 5 , k^l 2 ^l 3 ^U^l 5 
2 N ~ 4 , h = h , or, all permutations with 2 equal space indices 
2 N ~ 3 , h = h = h, or, all permutations with 3 equal space indices 
li = I2 and Is = 

or, all permutations 2 by 2 (n n n 

oN-2 7 7 7 7 nr [Ly.Z.l) 

Z , tl — <2 — '3 — '5? or 5 

all permutations with 4 equal space indices 

h = h = h and U = £5, or, 

all permutations with 2 or 3 equal space indices 
1 1 ; 7 1 ; 



= < 



2 N 1 , h = h = h = U = h- 



2) 



5N £ ni VJ 

Q 2 (h,l2,h,U,k) = / [[dr] I dr] I e 1 ^- 1 Vh Vh Vn+i 2 Vn+i 2 +i V2N+i 3 V2N+i 3 x 



1=1 

x V 3N+ i 4 i]sN+l A -l V4N+l 5 V4N+l 5 
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( 2 N 


-5 




-4 


2 N 


-3 



= < 



)AT-2 



h = h + 1 

£4 = + 1 and two other space indices are equal 
U = h + 1 and /1 = U and Z3 = £5, 

or, all permutations 2 by 2 (C.2.2) 

U = I2 + 1 and h = U = I5, or, 

all permutations with 3 equal space indices 

I2 = h — 1 = h and h = h = U. 



5N J2 *>' A Vj VJ 

/ JJ e 7 ' J=1 7^ 7^+^ ?7Ar+i 2 + l V2N+l 3 mN+h-i* 

J 7 1 



1=1 

x V3N+l 4 V3N+l 4 V4N+l 5 V4N+l 5 





h 


2 N ~\ 


h 


2 N ~\ 


h 



^2 + 1 and Zi = I4 and I3 = I5, or, 
< all permutations 2 by 2 (C.2.3) 

Z3 = h + 1 and h = I3 = U, or, 
all permutations with 3 equal space indices 
12^, h = l2 + ltmdl 1 = l s = U = l 5 . 



/ Yl driidfji e 1 '-^ 1 f) h r^-i fj N+ i 2 r] N+ i 2+1 f] 2N +i 3 r]2N+i 3 +iX 

J 7 1 



f nN-5 



)N-4 



X V3N+l 4 V3N+U-1 V4N+l 5 V4N+l 5 



h = h + 1 and I4 = I3 + 1, or, 
li = l 3 + 1 and I4 = I2 + 1 
h = ^2 + 1 and / 5 = / 3 + 1 and U = I5, or, 
Z2 = ^4 = Zi - 1 = Z3 + 1, or, 
h = h = h — 1 and Z 3 = Z 4 — 1, or, 

Zi = h + 1 and Z 4 = Z 5 = Z 2 + 1, or, (C.2.4) 
Z3 = Z5 = Zi — 1 and I2 = U — 1, or, 
Zi = Z3 = Z4 — 1 and I4 = I2 + 2 
Ar ~ 3 , Zi = Z3 = Z5 + 1 and Z2 = h and U = I2 + 2, or, 

Zi = Z5 H- 1 and Z2 = I5 and Z 3 = Z 5 — 1 and Z4 = Z5, or, 
Zi = Z3 = Z5 — 1 and I2 = h — 2 and Z4 = Z5. 
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5) 



5N 



VI A 



TT 

IJ VJ 



£5 (h,h,h,U,h)= J U drudrii e 1 ' 3 - 1 Vh+i Vn+i 2 Vn+i 2 V2N+i 3 V2N+i 3 x 

X V3N+l 4 V3N+l 4 -l VAN+l 5 V4N+l 5 



1=1 



( 2 N ~ 5 


h = 




u = 


2 N ~ 3 , 


h = 



= < 



■>N-2 



h + 1 

h + 1 and two other space indices are equal 
h + 1 and h = h = hi ° r ? 
all permutations with 3 equal space indices (C.2.5) 
1 4 = h + 1 and h = h and h = h, or, 
all permutations 2 by 2 
U = h + 1 and h = h = h = h — 1- 



6) 



{hi h, hi U, h) 



5iV 



r 52*11 A]\ 77.7 

J \\ driidfj! e 1 - 7 - 1 fj h r] h+1 fj N+h r] N +h-i V2N+i 3 V2n+i 3 +i>< 



1=1 



X V3N+U V3N+U-1 V4:N+l 5 V^N+ls 



2 N ~\ 



= < 



2 iv- 3) 

2 N-2 



h + 1 and U = h + 1, or, 
h + 1 and h = h + 1 
h + 1 and U = h = h + 1 
h = h + 1 and U = h + 1 
l 5 = h + 1 and h = h + 1 

1 



h 
h 
h 
h 
h 

h = U — 1 and h = h = h 
h = h = h — 1 and h = U 
h = h = h — 1 and h = U = h 



(C.2.6) 



7) 



{hi h-, hi hi h) 



5iV 



r 52 a Vj Vj 

/ II driidfji e 1 '^ 1 fj h rj h+1 fj N+ i 2 rjN+h+i V2N+i 3 V2N+i 3 -iX 

J 7 1 



7=1 



x V3N+U V3N+U-1 VAN+h V^N+h 
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/ oiV-5 



)N-4 



= < 



)iV-3 



h 


= h 


+ 1 and I4 = I2 + 1, or, 




h 


= h 


+ 1 and I3 = I2 + 1 




h 


= h 


+ 1 and I4 = I5 = I2 + 1, or, 




h 


= k 


= li + 1 and U = I2 + 1, or, 




h 


= k 


= I2 + 1 and U = l± + 1, or, 




h 


= k 


= li + 1 and h = h — 1? or 5 




h 


= h 


= Z 3 - 1 = h + 1, or h = l 3 = 


--h-l = l 2 + l 


h 


= h 


= 1$ and h = I5 — 1 and £3 = 


k + 1, or, 


h 


= h 


= I5 — 1 and £4 = /s and /2 = 


l 5 - 2, or, 


h 


= h 


= Z5 — 1 and ^3 = I5 and l\ = 


k-2. 



(C.2.7) 
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